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Abstract. Thermal history of the universe after big-bang nucleosynthesis (BBN) 
%jj I is well understood both theoretically and observationally, and recent cosmological 

r — ' observations also begin to reveal the inflationary dynamics. However, the epoch 

CNJ . between inflation and BBN is scarcely known. In this paper we show that the detection 

of the stochastic gravitational wave background around IHz provides useful information 
about thermal history well before BBN. In particular, the reheating temperature of 
^J I the universe may be determined by future space-based laser interferometer experiments 

QQ ' such as DECIGO and/or BBC if it is around 10^~^ GeV, depending on the tensor-to- 

f^ . scalar ratio r and dilution factor F. 
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1. Introduction 

Recent cosniological observations liave determined cosniological parameters with 
unprecedented accuracy [1]. From measurements of cosmic microwave background 
(CMB) anisotropy, the baryon content of the universe VLb is determined and it agrees with 
the prediction of big-bang nucleosynthesis (BBN). This confirms the standard thermal 
history of the universe at least for T < 0(1) MeV, where BBN begins. On the other 
hand, the primordial power spectrum of the density perturbation has been revealed to 
be nearly scale invariant. This has not only provided another evidence of inflation in 
the early Universe [2] but also given observational clues to study its details such as the 
shape of the scalar potential driving inflation [3J. 

Thermal history of the universe between inflationary era and BBN epoch, however, 
is least known observationally. In general, the universe is dominated by the coherent 
oscillation of the inflaton soon after the inflation ends, and finally the inflaton decays into 
standard model particles, which realizes the hot radiation-dominated universe. Thus 
understanding this "reheating" stage of the universe is very important. In a simple 
standard inflation scenario, the reheating effect is characterized by one parameter, the 
reheating temperature Tr |1], which is defined as 

where T^ is the decay rate of the inflaton, Mq = {8ttG)~^/'^ = 2.4 x 10^^ GeV is 
the reduced Planck scale and g^ denotes the relativistic effective degrees of freedom. 
However, Tr is hardly constrained from cosniological observations by now. Only the 
requirement is that Tr must be larger than a few MeV [5] in order that light elements 
must be created in a usual way. 

However, the reheating temperature Tr contains rich information from the 
viewpoints of particle physics. First, as is obvious from Eq. ([1]), Tr is determined 
by the inflaton decay rate, and it depends on the inflaton properties, such as its mass, 
potential and interaction strength with standard model particles. Thus determining 
Tr may have impacts on choosing realistic inflation models and inflaton candidates. 
Moreover, in supersymmetry (SUSY) [6], which is well-motivated physics beyond the 
standard model, theoretical upper bounds on the reheating temperature is imposed from 
the so-called gravitino problem [3, El El ^0\. The abundance of the gravitino, which is 
the superpartner of the graviton, is bounded from above in order not to destroy light 
elements through its decay processes for an unstable gravitino, or not to overdose the 
universe for a stable gravitino. This gives upper bounds on Tr, since gravitinos are 
produced efficiently in the reheating era and its abundance is proportional to Tr. 

Then, how can we probe such an early stage of the universe? As is well known, 
observed CMB photons come from the last scattering surface located at redshift 
z ~ 1100. The universe is opaque to photons beyond the last scattering surface, 
and hence we cannot directly look over the universe oi z > 1100 by observations of 
photons. Even the cosmic background neutrinos cannot be used as a probe of the early 
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universe with T > 1 MeV, since neutrinos couple with nucleons strongly in such a high 
temperature environment. 

The universe, however, is transparent to gravitational waves up to the Planck 
epoch in principle. This opens up the window for probing the very early universe, 
in particular the reheating epoch, with observations of gravitational waves. Along with 
many species of unwanted relics, inflation dilutes preexistent gravitational waves at the 
Planckian epoch. Instead quantum mechanical gravitational waves are produced in 
the inflationary era with an almost scale invariant power spectrum [llj. Although 
the amplitude of gravitational waves is constant in the super horizon regime, once 
a mode enters the horizon, it is reduced as the universe expands. Since the 
expansion rate depends on the equation of state of the universe, corresponding thermal 
history of the universe is imprinted in the gravitational wave spectrum at present. 
Although much work has been done which treats the spectrum and detection possibility 
of stochastic gravitational wave background of inflationary origin in the literature 

[la [131 [11 [ig [ini [HI [TBI [ini EDI [ai Ea Ea El ESI [^ 

did not take into account the modification of the gravitational wave spectrum around the 
reheating epoch and assumed nearly flat spectrum above the observationally interesting 
frequency ~ IHz, except for a few works where it was claimed that the equation of state 
of the early universe can be probed by looking at the spectrum of the gravitational wave 
background [21] (this kind of study was also done and extended in [2^). Recently, we 
have pointed out that taking into account the reheating effect on the gravitational wave 
spectrum is essential both for the detection itself and the purpose of probing thermal 
history of the universe between BBN and inflation, with an emphasis on its impacts on 
particle physics [3Tj . 

In this paper, we give more detailed analyses of reheating effects on the gravitational 
wave background and show that future mission concepts based on space laser 
interferometer, hke Japan's DECIGO [32] and NASA's Big Bang Observer (BBO), 
may have a chance to detect primordial gravitational wave background produced in the 
inflatonary era and determine the thermal history of the universe before BBN epoch, 
in particular the reheating temperature Tr. We also discuss that some "non-standard" 
cosmological scenarios, including late-time entropy production, can be probed. 

This paper is organized as follows. In Sec. [2]the spectrum of primordial gravitational 
wave backgrounds, taking into account the modification around the reheating epoch, is 
derived. Using this result, we discuss future prospects of the determination of the 
reheating temperature with space-based laser interferometers in Sec. [31 We discuss the 
case of late-time entropy production in Sec. [1 In Sec. [5l the impacts of determining the 
reheating temperature on particle physics models are summarized. Sec. [61 is devoted to 
our conclusions. 
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2. Primordial gravitational wave spectrum 

Primordial gravitational waves produced in the inflationary era have nearly scale- 
invariant spectrum and have effects on both large- and small-scale cosmological 
observations [12]. On large scales at observable frequeincies ~ 10~^^ Hz, the tensor 
metric perturbation generates B-mode polarization anisotropy of cosmic microwave 
background (CMB) for small multipole /, and such signal may be probed by next- 
generation ground-based or satellite observations of CMB (polarization) anisotropy 
[331 [m |35] . On small scale, gravitational waves may be detected by future space-based 
interferometer experiments with frequency around ~ IHz. Thus primordial gravitational 
waves give much information on wide range of the cosmological scales. In particular, 
direct detection of gravitational waves may reveal the state of the universe well before 
BBN, as we will see. 

Gravitational wave is described by the tensor perturbation on the metric, hij, which 
is deflned as 

cis^ = a\t) [-dr^ + {6ij + 2hij)dx'dx^] , (2) 

where hij is symmetric under the exchange of i and j, and satisfles the traceless and 
transverse condition, /i/ = 0, hij'-' = 0. Thus it has two physical degrees of freedom, 
which we denote as + and x. Conformal time r is deflned as dr = dt/a{t). Since tensor 
perturbation is gauge-invariant as it is, we do not need to care about gauge ambiguity 
as far as only the tensor mode is considered. The tensor perturbation hij is expanded 
in the Fourier space as 



A=+,x 



where the superscript A denotes each polarization degree, A = +/x, and e\, is the 



(A) 



polarization tensor, which satisfles e^^- e*-'*-^ ■* = S'^^ 

Primordial gravitational waves are produced in the inflationary epoch with an 
almost scale-invariant spectrum. The amplitude of the produced gravitational wave 
is proportional to the Hubble scale during inflation, -fTjnf . In terms of the dimensionless 
power-spectrum, it is given by [36] 



Af(fcr = 647rG f^ 



k I k 

l-2eln— + 2e(r/-e) (In — 



(4) 



counting two polarization states of the gravitational wave, where e and rj denote the 
slow-roll parameters during inflation deflned as e = Mq(V"'/1^)^/2 and rj = Mq{V" /V) 
with an inflaton potential V[(f)] and its flrst and second derivative with respective to the 
inflaton fleld 0, V and V". We take the pivot scale as A;* = 0.002 Mpc~^. Here slow-roll 
parameters should be evaluated at the epoch when CMB scale leaves the horizon. 
On the other hand, the spectrum of the curvature perturbation is given by 



2 n\ - 47rG /i/inf 
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+ (6e^ - 4er7 + 2r/^ - ( In 7^ 



k^' 



(5) 



where ^ = MqV V /V"^. Its normalization is already well determined observationally, 
as A|j ~ 2.0 X 10^^ on CMB scale [1]. Thus the slow-roll parameter e is related to 
tensor-to-scalar ratio r = A|/A^(A;*) through the relation 

r = 16e. (6) 

This shows that r is in general a small parameter and the tensor contribution is also 
expected to be small compared to the density perturbation. Notice that measuring r 
fixes the inflation scale as 

Vini = SM^Hlf ~ (3.2 X 10^*^ GeV)V. (7) 

Then let us derive the observable power spectrum of inflationary gravitational wave 
background. The evolution of the gravitational wave is described by the following 
equation, |j 

ht + 3Hht + ^ht = 0. (8) 

The amplitude of the gravitational wave with comoving wave number k remains constant 
when the mode lies outside the horizon. However, once it entered the horizon, its 
amplitude begins to damp. For the mode which enters the horizon in the matter- 
dominated regime, the solution is written as 

I4(r)=ht^"[^), (9) 

where je denotes the i-th spherical Bessel function. In general, the solution in a power 
law background a{t) oc f^ can be written in the form 

/ifc(r) oc a(t)^ J 3p-i (kr), (10) 

2(l-p) 

with Bessel function Jn{x). Another damping factor comes from the fact that the 
relativistic degrees of freedom g^, does not remain constant and expansion rate is modifled 
from simple power law a{t) oc T^^ in the early universe. This effect gives the damping 
factor [281 ES] 



I \ (T M ' ^^^) 

9*0 / \9*s{'L in) / 

on the power spectrum of the gravitational wave, where Tin denotes the temperature at 
which the corresponding mode enters the horizon, given by 

' 9*s{Ti ^^ '^''^ 
106.75 7 VlO^^Mpc 



Tin ^ 5.8 X 10« GeV ( '^^^ ] ( ,^,,:,.^„„, ) . (12) 



I Here we neglect the anisotropic stress. Including an anisotropic stress modifies the gravitational 
wave spectrum for the present frequency < 10^^ Hz due to the neutrino free streaming effect [571I29J . 
but our concern is the frequency around ^ 1 Hz. Hence the following arguments are not modified. 
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Moreover, recent cosmological observations have revealed that the expansion of the 
present universe is accelerating due to the dominance of unknown energy density, called 
dark energy (here we assume cosmological constant for simplicity). This also gives 
suppression factor ~ {flrn/^AY on the power spectrum |16j. As a result, the present 
gravitational wave spectrum per log frequency interval is written in the form 



2 



figw(/) = Y^^^lik), (13) 



where 



Alik) = A 



(P)fU\'i 




x,^)T^ixn), (14) 



with a bar denoting the average over many periods. Here, Ti(a;eq) denotes the transfer 
function, which connects the gravitational wave spectrum of the mode which enters the 
horizon before and after matter-radiation equality, t = tgq. It is calculated as [16] 

T2(a;eq) = [1 + 1.57xeq + 3.42xy , (15) 

where Xeq = k/k^q and keq = a{tcq)H{teq) = 7.1 X 10~^f2m/i^ Mpc~^|§| 

On the other hand, T2{xr) connects the mode which enters the horizon after and 
before the reheating ends. Because before the inflaton decays the universe is dominated 
by the coherent oscillation of the inflaton, the spectrum of the gravitational wave changes 
for the mode which enters the horizon at the infiaton-dominated epoch [k > fc^), where 

kn ^ 1.7 X 10- Mpc- f'igi) ''' (-^) . (16) 



^106.75; VlO^GeV^ 
In terms of the frequency, this corresponds to 

;,. 0.026 H.(^)"(,^). (17) 

which is close to the most sensitive frequency band of DECIGO and BBO for Tr ~ 
10® GeV. The transfer function T2{xr) is obtained by solving simultaneously Eq. ([8]) 
and the following Friedmann equations taking the decay of inflaton into acount, 

P0 + 3Hp^ = -T^p^, (18) 

Pr + AHpr = T^p^, (19) 



H' = ^ip, + Pr), (20) 

where p^ denotes the energy density of the inflaton coherent oscillation. We flnd that 
T2{xji) is well approximated by 

T^ixR) = [1 - 0.32XR + 0.99xy '' , (21) 

§ CoefBcients in the right hand side of (jisp differ from those of Ref. [16] reflecting the difference of 
the definition of fcoq. 
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Figure 1. A schematic picture of evolution of the Hubble horizon scale H^^ 
and physical wave length of some modes. (/).D., R.D., M.D. and A.D. denote the 
inflaton oscillation dominated era, radiation dominated era, matter dominated era, 
and cosmological constant dominated era, respectively. 



where xr = k/k^. Thus Qg^{f) behaves as oc f~'^{f^) for the mode which enters in 
the horizon in the matter (radiation) dominated regime. In Fig. [1] we show a schematic 
picture for the evolution of the Hubble horizon scale H~^ by solid line, and physical 
wave length of various modes, k > kji, fceq < k < kji and k < fcgq by dotted lines. 
Since the evolution of inflationary gravitational wave is sensitive to the expansion of the 
universe, H~^, Fig. [T]shows that characteristic feature of power spectrum can be probed 
for the mode around / ~ f^. 

Before discussing the detection possibility, we briefly mention the effect of running 
spectral index. The running spectral index depends on the slow-roll parameters e and 
rj. Although e is related to r through the relation ([6]), we are free to choose r] as long 
as the scalar spectral index ns{— 1 — 6e + 2ri) lies in the favored range from recent 
WMAP results. In Fig. [2] the resulting spectrum flg^{f) at / = 0.1 Hz is shown for 
T] = 0.01, 0, —0.01 from upper to lower. Here sufficiently large reheating temperature 
Tr ^ 10^ GeV is assumed. As is clearly seen, for r > 0.1 the gravitational wave 
amplitude at / = 0.1 Hz decreases since the spectral tilt and its running becomes large. 
The difference of rj yields slight change on the gravitational wave amplitude for r < 0.1 
(c.f. the current contraint on r is r < 0.20 (95% CL.)). Hereafter we set t] = for 
simplicity. 

In Fig. [3l we show the resulting gravitational wave spectrum for Tr = 10^ 
and 10^ GeV with r = 0.1 and 0.001. Together with the predicted gravitational 
wave spectrum, sensitivities of the future space-based laser interferometer experiments, 
DECIGO, correlated analysis of DECIGO, ultimate-DECIGO (single) and correlated 
analysis of ultimate-DECIGO [S^ are also presented. It can be seen that theoretical 
predictions are well above the sensitivity line of DECIGO (correlated) and ultimate- 
DECIGO for large enough Tr and r, while the direct detection is undesirable for 
Tr < 10^ GeV. 
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3. Prospects for the determination of Tr with future space-based laser 
interferometer experiments 

The spectrum of the primordial gravitational wave background generated during 
inflation crucially depends on the reheating temperature Tr after inflation, as can be 
seen from Fig. [3l Conversely, this fact opens up the possibility that future experiments 
devoted to detect gravitational wave background will probe the reheating stage of the 
universe. 
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Figure 2. J7gw(/) at / = O.lHz for 77 = 0.01, 0, —0.01 from upper to lower. 
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Figure 3. Primordial gravitational wave spectrum for Tr = 10^ GeV and Tfi = 
10^ GeV are shown by thin and thick lines for r = 0.1 and 0.001. Also shown are 
expected sensitivity of DECIGO (green dashed), correlated analysis of DECIGO (blue 
dot-dashed), ultimate-DECIGO (purple dashed) and correlated analysis of ultimate- 
DECIGO (red dotted), from upper to lower. 
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The important parameters that determine the primordial gravitational wave 
spectrum are the tensor-to-scalar ratio r and the reheating temperature T^. Tensor- 
to-scalar ratio determines the overall normalization of the spectrum, and Tr fixes the 
frequency above which the spectrum is significantly suppressed. The important point is 
that if the bending point of the gravitational wave spectrum determined by Tr around 
the frequency given by flT6|) lies above the sensitivity of detectors, Tr can be determined 
by observations of gravitational waves. 

The non-zero value of tensor-to-scalar ratio r will be probed with measurements 
of the B-mode CMB polarization [31]. Since the B-mode is generated by the tensor 
perturbation only for large angular scale, i < 100, detection of the B-mode polarization 
indirectly confirms the existence of the gravitational wave background. Notice that this 
effect is seen on very large-scale anisotropy with wavelength comparable to the present 
horizon scale, so the detection of the B-mode polarization is somewhat complementary to 
the direct detection of the gravitational wave on small scale with wavelength of the order 
of the Earth radius. The Planck satellite |39j, scheduled to be launched in October 2008, 
will measure r up to ~ 0.1. On-going ground-based or balloon experiments devoted to 
detect CMB polarization, such as Q/U Imaging ExperimenT (QUIET) |10] and Clover 
[H] will detect r up to ~ 0.01. Ultimately future space mission will dedicate to detect 
the primary B-mode signal up to r > 10~^ [34]. According to recent works [ISIHS], the 
observed value of the scalar spectral index, Ug ~ 0.961 ± 0.017, implies somewhat large 
tensor contribution, r > 10~^ for various inflation models. Thus it seems to be plausible 
that non-zero r will be conflrmed by CMB polarization experiments. 

Once r is measured from those experiments, we may have a chance to detect signals 
of gravitational wave background of inflationary origin. Planned Laser Interferometer 
Space Antenna (LISA) [Hj cannot reach the required sensitivity even the largest possible 
value of r allowed by WMAP5 and high enough Tr are assumed. However, future 
mission concepts like DECIGO and BBO are likely to detect them or put a meaningful 
constraint on Tr. Therefore our main interest is to demonstrate the observable range of 
reheating temperature through such future experiments for various values of r. Before 
discussing results, however, it should be noticed that stochastic noise coming from 
some astrophysical processes must be taken into account for the purpose of detecting 
primordial gravitational wave background. In particular, gravitational waves from 
white-dwarf binaries are considered to completely hide the primordial ones for the 
frequency / < O.lHz [15]. But for the frequency range 0.1-lOHz, where the DECIGO 
and BBO are most sensitive, foregrounds from astrophysical objects are separable. [||| 

In Fig. m we show parameter region on r-T^ plane which will be probed with 
the correlated analysis of DECIGO, ultimate-DECIGO (single) and ultimate-DECIGO 
(correlation) respectively. We have cut the sensitivity of these projects below / ^ 0.1 Hz, 

II Gravitational waves from the collapses of Population III stars may be a dominant contribution 
around the deci-Heltz band [46 . However, it crucially depends on the early star formation rate, which 
is very uncertain now, and such a signal would be separable if we adopt reliable abundance of Population 

III stars and take the duty cycle and their angular distribution into account. 
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taking into account the stochastic noise from white-dwarf binaries. The gravitational 
wave background can be detected in the hght blue shaded region. Furthermore the dark 
blue shaded region shows the parameter region where the value of Tpi can be determined 
with signal-to-noise ratio 5. It is seen that for 10~^ ^ '^ ^ 1; direct detection of 
gravitational wave background can determine the reheating temperature Tr, if it lies in 
the range Tr ~ 10^-10^ GeV. 

If higher Tr is realized in nature, the inflationary gravitational wave background 
can be detected by DECIGO or BBO, but the value of Tr remains undetermined. In 
this case, the lower bound on Tr can be read off from these figures. For example if 
r = 0.1, the detection of gravitational wave background means Tr > 10^ GeV, which 
provides useful constraints on some particle physics models, as we will see in Sec. |5l 

4. Late-time entropy production 

So far, we have assumed that there were no late-time entropy production processes 
after the completion of reheating after infiation. However, this may be too simplified 
assumption. Let us consider the case some scalar field x other than the infiaton 
dominates the universe after the infiaton decays and x eventually decays releasing 
huge entropy. Examples are Polonyi [17] or moduli field [HI HO], scalar partner of 
the axion [50], or others. Such late-decaying particles are interesting since they dilute 
cosmologically harmful gravitinos. We denote the dilution factor as F, defined by 
^ s{T^)a^{T^) _Tr ( p^^ 



s{Tn)a%TR) T^ Kp^Jt^' 

where T^ is the decay temperature of x, which must be larger than a few MeV, and 
p^ denotes the energy density of the x-field coherent oscillation. [^D For a simple model 
where the x-field has an initial amplitude Xi ^"^^ begins to oscillate when the Hubble 
parameter becomes equal to the mass of the x, it is estimated as 

r. Tn xl 



T,3M^- 



(23) 



The abundance of all dangerous cosmological relics produced in the reheating era after 
infiation, such as the gravitino discussed below, are diluted by this factor. 

Importantly, such non-standard cosmological evolution scenarios are imprinted in 
the present gravitational wave spectrum. In the presence of such a late-decaying particle, 
the additional x-matter dominated era suppresses the gravitational wave amplitude for 
the frequency which reentered the horizon during or before x begins to dominate the 
universe. (Fig. [5] shows a schematic picture.) The spectrum now becomes 

n^if, F) = n,M) X T^ix^)T^{x^R), (24) 

where flg^{f) is given by Eq. (TT3l) with k^ replaced by 

kniF) = kRF-^'\ (25) 

^ Here we assume that x begins to oscillate during the infiaton oscillation dominated phase. 
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Figure 4. In the outer light shaded region the gravitational wave background can 
be detected, and the inner blue shaded region shows the region where Tr can be 
determined with signal-to-noise ratio 5 by correlated analysis of DECIGO, ultimate- 
DECIGO (single) and ultimate-DECIGO (correlation) from upper to lower. 
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Figure 5. A schematic picture of evolution of the Hubble horizon scale _ff~^ and 
physical wave length in the presence of late-time entropy production from x- X-D- 
represents x-dominated era. 



Here Xy{= k/k^ corresponds to the wavenumber which enters the horizon at the decay 
of X, 
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(26) 



and x^r{= k/k^fi) corresponds to the epoch when ^-domination begins, given by 

k^n{F) = k^F^/\ (27) 

We can see that for the mode k^n < k < kji, which corresponds to the mode which 
reenters the horizon in the radiation dominated era before the x- domination, the energy 
density of the gravitational waves is suppressed by the factor ~ (A;^/A;^^)^ = F~^/^ [2T]. 
On the other hand, there are no effects on large scale with the mode k < k^. 

Thus the gravitational wave spectrum in the presence of late-time entropy 
production is completely characterized by two additional parameters, the dilution factor 
F and the decay temperature of x, T^- However, in many cases T^ is expected to be 
very small, say, ~ 0(1) MeV-O(l) GeV, and hence it does not affect the gravitational 
wave amplitude at the frequency relevant for the direct detection. Thus hereafter we 
mainly focus on the effect of varying F. 

Note that non-negligible F affects only the overall amplitude of the gravitational 
wave background for the mode k^R < k < kn, and hence there is a degeneracy between 
F and tensor-to-scalar ratio r when concerning the direct detection around O.lHz. 
However, the tensor-to-scalar ratio should be determined by cosmologically large scale 
CMB B-mode anisotropy. Thus if future CMB experiments measure r, there does not 
remain an ambiguity coming from the degeneracy between r and F. In other words, 
if the result of direct detection deviates from the expected signal from the large scale 
measurement of r, there must be an entropy production process in the early universe. 
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All of these features are seen in Fig. El where the gravitational wave spectrum 
with F = 10^ and 10^ are shown. Here we have fixed r = 0.1, Tr = 10^ GeV and 
Tj(, = 1 GeV. This figure shows how the gravitational wave spectrum is affected by the 
late-time entropy production. 
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Figure 6. Gravitational wave spectrum for the dilution factor F ~ lO'^ and 10"'. Here 
we have fixed r = 0.1, Tr = 10^ GeV and T^ = 1 GeV. 

If 0.1 Hz< kji{F) < 10 Hz, both F and Tr can be determined from the shape of 
the gravitational wave spectrum. In Fig. [6] we show future sensitivity to determine the 
dilution factor F and the reheating temperature Tr with fixed tensor-to-scalar ratio r, 
which can be measured by CMB polarization experiments. In Fig. [7] we show parameter 
region in F—T^ plane where both parameters can be determined by correlated analysis of 
DECIGO, ultimate DECIGO (single), and ultimate DECIGO (correlation) respectively 
by the dark shaded region. The light shaded region shows the parameters where the 
gravitational wave background is detected, but the value of Tr remains undetermined. 
Here we have fixed r = 0.1 and T^ = 1 GeV. (The precise value of T^ does not matter 
as long as T^ < 10^ GeV.) On the other hand if fc^(F) < 0.1 Hz, we can measure only 
the ratio T^/F. However, this ratio contains sufficient information to determine the 
gravitino abundance including the dilution effect. Finally if /cr(F) > 10 Hz, only a 
lower bound on Tr is obtained as Tr > 2 x lO^F^/^ GeV. Even this case is useful for 
constraining the gravitino mass, as we will see. 

5. Implications on particle physics 



We have shown that the reheating temperature of the universe Tr may be determined 
from future space-based laser interferometer experiments. In this section we discuss 
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Figure 7. In the outer light shaded region the gravitational wave background can 
be detected, and the inner blue shaded region shows the region where Tji can be 
determined with signal-to-noise ratio 5 by correlated analysis of DECIGO, ultimate- 
DECIGO (single) and ultimate-DECIGO (correlation) from upper to lower. Here 
r — 0.1 is assumed. 
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some implications of determining Tr on particle physics, in particular SUSY modelilJ 
and baryogenesis mechanisms. 

5.1. Supersymmetric models 

5.1.1. Thermally produced gravitinos In SUSY, there exists the superpartner of the 
graviton, the gravitino. The gravitino mass ranges from 0(1) eV to O(IOO) TeV 
depending on SUSY breaking models, while other SUSY particles have the mass of 
0(1) TeV. Gravitinos are produced in the early universe through scattering of particles 
in thermal bath, and the resultant abundance of the gravitino is proportional to Tji 
[81EI1E2], 



where Y3/2 = 'n-3/2/s is the gravitino number-to-entropy ratio, 7713/2 denotes the gravitino 
mass and mg(~ 0(1) TeV) is the mass of the gluino, fermionic superpartner of the gluon. 
If the gravitino is rather heavy {ms/2 > 1 TeV) and unstable, it eventually decays with 
lifetime typically longer than 1 sec producing high-energy photons and hadrons. Those 
decay processes destroy or overproduce light elements such as ^He, D, ^Li and ^Li. This 
constrains the reheating temperature as Tr < 10^^^ GeV depending on the gravitino 
mass and its hadronic branching ratio [8j. On the other hand, if the gravitino is the 
lightest SUSY particle and stable due to the i?-parity conservation, it contributes to 
total matter density of the universe [10] . This leads to the constraint 

.,<,.,0»Gev(i^)(^)-V. (., 

for m3/2 ~ 10-^-10 GeV. 

For example, let us consider the situation where Tr is revealed to be larger than 
~ 10^ GeV and negligible dilution factor {F ~ 1) is confirmed with future space- 
based laser interferometer experiments. For the unstable gravitino with 7/13/2 ^ 1 TeV, 
as is often the case with gravity-mediated SUSY breaking models, its abundance is 
constrained as F3/2 < 10^^^ |8j, and hence Tr should be less than around ~ 10® GeV. 
Thus SUSY breaking models which predict the gravitino mass of 0(1) TeV will be 
excluded, even if the gravitino mass might not be determined by accelerator experiments. 
Also for the light gravitino scenario, the following constraint is obtained, 

my. > 1 GeV (^] (jj^) i^y. (30) 



^Fj \W GeV J Vl TeV. 
Thus the gravitino mass with 77/3/2 ^ 1 GeV, which can be realized in some classes of 
gauge-mediated SUSY breaking models [53], is excluded since otherwise the gravitino 
abundance exceeds the present dark matter abundance. 

+ SUSY theories contain nearly twice particle species compared with the standard model and the 
relativistic effective degrees of freedom at high temperature is doubled, g*{T > 1 TeV) = 228.75. This 
leads to slight suppression on the gravitational wave spectrum, but does not much affect the results of 
SecEl 
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On the other hand, if the gravitino is stable and its mass is determined from 
accelerator experiments, Tr is the only parameter which determines the present gravitino 
abundance. If the observationally inferred value of Tr saturates the upper bound ( l29l) . 
the gravitino abundance is appropriate for the dark matter of the universeo Note that 
if the dark matter consists of the gravitino, it is impossible to detect it directly because 
of its too weak interaction strength. Thus Tr has an important information to probe 
whether the gravitino truly takes a roll of the dark matter or not. 

5.1.2. Non-thermally produced gravitinos Recently it is pointed out that non-thermal 
production process of gravitinos from the inflaton decay gives a significant amount of 
gravitino abundance [56], 

where m^ and (0) denote the mass and VEV of the inflaton. Since non-thermal 
contribution is proportional to T^^, this gives lower bound on Tr for fixed m^i2. 
Furthermore the inflaton always decay into MSSM sector appearing in the superpotential 
in supergravity, which provides lower limit on the reheating temperature as [57] 

T >inTVl I /^ 228.75 y/y m^ ^^ f (0) ^ 

where yt is the top Yukawa coupling. 

In Fig. [8] allowed parameter region in m^i2-TR is shown for the case of stable 
gravitino for m^ = 10^'^ GeV and (0) = 10^^ GeV. The upper left and lower right 
regions are excluded from thermal fl28|l and nonthermal production fl3Tl) . respectively. 
Also spontaneous decay of the inflaton gives lower bound on Tr ( !32l) as denoted by the 
dotted line. It is seen that 10^ GeV ^ Tr < 10^ GeV is favored, and this range is 
interesting from the viewpoint of direct detection of gravitational waves. 

5.2. Baryogenesis mechanism 

Another important issue which is deeply related to the reheating temperature is 
baryogenesis. As is well known, almost all the ordinary matter in the universe consists of 
baryon, not anti-baryon. It is a long-standing mystery how the sizable amount of baryon 
asymmetry is generated after inflation. We must rely on some new physics which involves 
extra baryon number violation, CP violation and non-equilibrium process in order to 
create matter-anti-matter asymmetry. One of the most popular mechanisms to create a 
correct amount of baryon asymmetry is thermal leptogenesis scenario using right-handed 
neutrino [58]. In order for this mechanism to work well, Tr > 10^ GeV is required |59j . 
Right-handed (s)neutrinos can also be produced non-thermally via right-handed 
sneutrino condensation [60l [61] or inflaton decay [621 E3] , and Tr > 10^ GeV is required 



* Gravitinos produced by the decay of the next-to-hghtest SUSY particle may also be the dark matter 
[54l [55] , but BBN constraints almost exclude such a possibility. 
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in order to create observed amount of baryon asymmetry in these scenarios. Thus non- 
thermal leptogenesis scenarios can also be favored or disfavored from future space-based 
gravitational wave detectors. 

Besides those leptogenesis scenarios, there are many other baryogenesis mechanisms 
which we do not list here [6l], and many of them predicts baryon asymmetry 
proportional to Tr. Thus determining Tr has important implications on baryogenesis 
mechanisms. 



6. Conclusions and discussion 

In this paper we have shown that direct detection of stochastic gravitational wave 
background of inflationary origin carries rich information on the early universe. 
In particular, the reheating temperature of the universe after inflation Tr can be 
determined with future space-based laser interferometer experiments, DECIGO and 
BBO. In case that non-zero tensor-to-scalar ratio is confirmed by CMB polarization 
observations, DECIGO or BBO will detect signals of gravitational wave background if 
Tfj > 10^ GeV. Moreover, they can not only detect gravitational waves, but determine 
Tr if Tr < 10^ GeV, through the Tr dependence of gravitational wave spectrum. 
Thus combined analysis of direct detection of gravitational wave background and CMB 
anisotropy measurements provide consistency check of inflation models. It also restricts 
realistic model among many and many inflation models, since both the tensor-to-scalar 
ratio and the reheating temperature crucially depend on the properties of the inflaton - 
its mass, potential, interaction strength, etc. 
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We have also discussed implications of determination of the reheating temperature 
on particle physics. The detection of gravitational wave background clearly goes beyond 
just a probe of thermal history before BBN. In the environment of the very early universe 
with extremely high temperature, currently undiscovered particles predicted by some 
physics beyond the standard model may be efficiently produced. Thus probing this 
epoch is directly connected to particle physics. For example, some SUSY breaking 
models will be excluded or severely constrained from the so-called gravitino problem, if 
future observations determine Tr. Also some baryogenesis scenarios, including thermal 
leptogenesis scenario, may be excluded or disfavored. In SUSY axion models, constraints 
on Tr may become much more stringent [65l |66] . 

As a final remark, gravitational waves may also be generated from some other 
cosmological processes, such as preheating after infiation j67l[6Hl[69] and first order phase 
transition followed by subsequent bubble collisions [70l[7l]. Although the amplitude and 
typical frequency from these contributions are highly mo del- dependent, they might give 
complementary information on infiation models or cosmological evolution scenario after 
infiation, if detected independently of primordial gravitational waves analyzed in this 
paper. Moreover, recently it is pointed out that decay of domain walls associated with 
gaugino condensation may produce a significant amount of gravitational waves and can 
be used as a probe of the gravitino mass 
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